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Abstract 



Attention is drawn to the observation that in many correlated systems (e.g. 
\ 3 He, heavy fermion systems and Hubbard models) the specific heat curves, 

o, 

when plotted for different values of some thermodynamic variable (e.g. pres- 
et 



sure, magnetic field, interaction), cross almost precisely at one or two charac- 
teristic temperatures. A quantitative explanation of this phenomenon, based 



Sh , on the temperature dependence of the associated generalized susceptibilities, 

O 

is presented. 

> ■ 

^ ■ PACS: G7.55.Cx. 7J..27.+a.7L2S. + rl 



In 1959 Brewer et al. |l| noticed that the specific heat curves C(T,P) of liquid 3 He 
measured at different pressures P all intersect at a temperature T + ~ 0.16i^, and that 
C(T,P) increases with increasing pressure below T + . Greywall's high-precision measure- 
ments of the specific heat at constant volume V provided striking evidence for the sharpness 
of the crossing point at T + ~ 160mK (see Fig. la) @. At this temperature the specific heat 
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is obviously independent of volume or pressure. It was unclear, however, whether special 
significance should be attached to this finding 0,0. Recently, Georges and Krauth [f| ob- 
served the same conspicuous crossing phenomenon in quite a different system, namely in the 
paramagnetic phase of the Hubbard model, the simplest model of correlated electrons, in 
infinite dimensions. For small to intermediate values of the local interaction U the specific 
heat curves C(T,U), calculated by iterated perturbation theory, were found to intersect 
almost at the same temperature T + ~ 0.59t*, where t* is the scaled hopping amplitude of 
the electrons (Fig. 2). Clearly the existence of these peculiar points of intersection calls for 
an explanation. 

In this Letter we illustrate that crossing points such as the ones described above can 
actually be observed in the specific heat of many correlated systems. Furthermore, we 
explain the origin of this phenomenon. 

To be able to discuss the problem in a sufficiently general frame work we define a gen- 
eral free energy $(T, X), where X can be any thermodynamic variable, e.g. pressure (P), 
magnetic field (P), on-site interaction (£/). The conjugate variable associated with X is 
£(T,X) = -d$(T,X)/dX. Pairs of variables (X,£) are, for example, (P, — V), (P, M) and 
(U, —D), with M as the magnetic moment, and D as the number of doubly occupied sites in 
Hubbard models. With the entropy S(T,X) = — <9$/<9T one obtains the Maxwell relation 

dS(T,X) _ d£(T,X) 

dx dr ' u 

A search of the literature reveals that there exists quite a number of systems, both in 
theory and experiment, where the specific heat curves C(T,X) = TdS(T,X)/dT versus T 
when plotted for different, not too large values of X intersect at one or even two well-defined, 
non-zero temperatures. Apart from normal-liquid 3 He it can be observed in heavy fermion 
systems with and without Fermi liquid behavior, for example in CeAl 3 || (Fig. 3a) and 
UBei3 upon change of P, in UPt3_ x Pda; || and CePtaSii-zGe^ as x is varied, and 
in CeCug-xAu^ (x = 0.3,0.5) when either P [ID] or B [|H{|(Fig. 3b) is varied. It is also 



found in Euo.5Sro.5As3, a semi-metal with competing interactions, upon change of B JT2 



In particular, all theoretical models of Fermi- and Luttinger liquids investigated beyond the 
low-temperature regime show this features: the one-dimensional (d = 1) Hubbard model in 



a magnetic field ||13||, the 1/r-Hubbard in d = 1 in the metallic phase when the interaction 



U is changed |TJ]] , and the Hubbard model in d = oo discussed above Q . 

To explain the origin of the crossing points we separate the problem into two questions: 
(i) Why do specific heat curves cross at all?, and (ii) How wide is the region where the 
curves cross? Turning to the first question, we note that any crossing of specific heat curves 
C (T, X) implies 



dC(T, X) 



d 2 £(T,X) 



T+(X) oT 2 



= 0. (2) 

T+(X) 



OX 

Thus crossing occurs where £(T,X) versus T has a turning point. In general the crossing 
temperature T + (X) still depends on X. Only if T + is independent of X for some range 
of X- values do the curves intersect at one point. Crossing of specific heat curves may be 
inferred from a sum rule for the change of the entropy S(T,X) with respect to X in the 
limit T — > oo 

! dS(T,X) X f °°dT' dC(T',X) 

1. Lattice models (X = U): Eq. (g) implies r\ v = for any kind of Hubbard model since 
S approaches a constant for T — > oo. (At high temperatures, T ^> U,C(T,U) oc U/T, i.e. 
dC/dll > 0). Hence dC/dU must have positive and negative contributions, i. e. the specific 
heat curves must cross at one or more temperatures, for the integral to vanish identically. We 
note that this is a genuine correlation effect originating from the existence of U 2 and higher 
terms in an expansion of the internal energy E(T, U), and hence of C(T, U) = dE/dT, in 
powers of U. 

2. Continuum systems (X = P): Here eq. (Jfy implies rj P = —1 since S approaches the 
ideal-gas value for T — > oo. Apparently dC/dP < at most (especially high) temperatures 
T5fl. Now, in the Fermi liquid phase of 3 He JDS] and in the paramagnetic phase of the 



Hubbard model at low temperatures [fH,|l3| the entropy is known to increase with X(= P, U), 



i.e. BC/dX > 0. Hence, at sufficiently low T when S(T,X) = C(T,X) = 7T, one has 
d^/dX > 0. This may be attributed to the excitation of low-energy (spin) degrees of freedom 
in the correlated system ||16|| . Eq. (3) then implies that the specific heat curves of liquid 3 He 
cross once (Fig. la) and in the Hubbard model twice (Fig. 2). These two systems only consist 
of a single species of particles and there exists one characteristic low-temperature scale - the 
crossing temperature T + (X). (In fact, in 3 He T + practically coincides with the temperature 
above which Fermi liquid theory breaks down fl7|). By contrast, heavy fermion systems 
are basically two-component systems consisting of conduction and localized electrons which 
may hybridize. The strength of the hybridization is determined by an amplitude Vh y b{P) 
which increases with pressure. By hybridizing, the electrons may gain an energy fc^T^, 
where T K (P) oc exp(— const /Vfi b ) is a different ("Kondo") low-energy scale. Below T K the 
linear specific heat coefficient is given by j(P) oc where now d^/dP < 0, and hence 

the specific heat decreases with pressure. This implies that the specific heat curves will cross 
twice: at T + , below which the low-energy (spin) excitations lead to dC/dP > 0, and again 
at T' + < Tk below which dC/dP < 0. This is precisely what is seen in several heavy fermion 
systems, e.g. in CeAl 3 [§ (T' + ~ 5K, T + ~ 17K; see Fig. 3a) and UBe 13 [§ [T' + ~ 2.5K, 
T + ^9K) @. 

We now turn to the question concerning the width of the crossing region. A plot of 
C(T, X) versus X (see e.g. Fig. lb for 3 He) shows that C(T,X) depends only weakly on X 
even for large values of X. Hence we expand C(T, X) in an (asymptotic) series in X — X , 
with X Q chosen at convenience, 



C(T,X)~C(T,X ) + (X-X )T 



dT 



x 



\x-X,fT d * 



+ ... (4) 

X 



2 v uy q T 2 Q X 

where we used eq. (2). At T+(X ) eq. (|) implies C(T+,X) ~ C(T + ,X )[l + W Xo (X)] 



where the (relative) width of the crossing region, |Hx (X)| = |A^ (X) + A^(X) + . . . |, is 



determined by the numbers 
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A ^ (X) -(n + l)!C(T + ,X ) df^ X (T ' X ' 



2 

(5) 



with x {n) (T,X) = d n £/dX n . For |A^(X)| < 1 the C(T,X) curves will intersect at a 
well-defined point. The width is seen to be determined by the curvature (with respect to T) 
of the linear (n = 1) and non-linear (n > 1) susceptibilities x^ n \T,X) at T + and Xq. There 
are two particularly relevant cases in which the A^ n ^ are small: 

(i) Weak T-dependence of x ■ F° r X = P, £ = —V the susceptibility x'' 1 '* — —dV/dP = 
K^V is essentially the isothermal compressibility of the system. For 3 He and heavy-electron 
liquids we can estimate Ap^ (P) from the linear-T behavior of the specific heat which, using 
eq. (|I|), implies x (n) ( T i P ) - X in) (0,P) + \^ n+1 \P)T 2 , with ^ n \X) = d n ^/dX n . Thus 
we obtain Ap^(P) ~ [(1 - P/P )" +1 /(n + l)l] [P n+1 7 {n+1 )(P )/7(P )] . For 3 He (with 
P = 15bar,P 2 7 ( 2 )(P )/7(Po) ~ 5 x 1(T 2 ) we find \W(P)\ < 0.03 for < P < 30 bar, i.e. 
on the scale of Fig. la the crossing region is indeed essentially confined to a point. Hence 
it is the very weak temperature dependence of —dV/ dP of liquid 3 He at low temperatures, 
i.e. the almost linear pressure dependence of the 7-factor of the specific heat, which is the 
origin of the sharp (but not exact) crossing of specific heat curves at T + ~ 160mK. Similarly, 
for the crossing at T' + ~ 5K in CeAl 3 [g] (with P = 4.8 kbar, P o y 2 )(P )/7(P ) ~ 0.4) we 
obtain |W(P)| 0.2 for P between 0.4 — 8.2 kbar, which again agrees with the data (Fig. 
3a). 

(ii) Linear X-dependence of x^'- If £ is a linear function of X, i.e. £(T,X) = x^{T)X, as 
in linear-response theory, the crossing condition, eq.(0), takes the form d 2 x(T + )/dT 2 = 0, 
where T + = T + (X = 0). This implies that Aq^ vanishes identically for n > 1. In this 
case all specific heat curves intersect exactly at one point. The width |Wo(X)| becomes 
finite only through non- linear terms in £(T,X) = x {1) (T)X + ±x {3) (T)X 3 + . . ., where 
■yA n )(T) = x^ n \T,X = 0). The lowest-order contribution to the width is A^(X) which 
may be quite small. This explains why for not too large magnetic fields the specific heat 
curves C(T,B) cross at a well-defined temperature both in the d — 1 Hubbard model 
at U — const |13|, as well as in CeC^-zAu^ ||TT| , [T9|1 (Fig. 3b). The same arguments 



apply to C(T, U) of the paramagnetic phase of Hubbard models p]Jl4}| where we now choose 
£ = -D(T, C/) = | — -D(T, Z7) at half filling such that D(T, 0) = 0. To a good approximation 
D(T, U) is linear in U for not too large U at all temperatures P0| . We find |W(i7)| < 0.05 
for U < 2.5t* in the d — oo Hubbard model ]2T| . 

In Hubbard models the intersection of C(T,U)/kB = f(T/t,U/t) curves is sharp only 
at high temperatures. At low temperatures the generation of low-energy excitations leads 
to a renormalized energy scale t — > t e jj^i t. Hence a perturbation expansion of E(T, U) or 
C (T, [/) to second order in U will be valid only for a small range of [/-values, implying a 
wide crossing region, |W(Z7)| ~ 1, at low temperatures. 

In summary, we showed that the remarkable crossing of specific heat curves C(T,X) 
vs. T for different thermodynamic variables X, first observed in 3 He @ and the Hubbard 
model Q , is not accidental but can be found in many correlated systems. The width of the 
crossing region is found to give explicit information about the temperature dependence of 
the generalized susceptibilities associated with X. A related observation, that for Hubbard 
models the value of the specific heat at the crossing point is almost universal, will be discussed 



elsewhere 22 
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FIG. 2. Specific heat C(T,U) of the paramagnetic phase of the Hubbard model in d = oo 
dimensions calculated by iterated perturbation theory ||. 

FIG. 3. Specific heat (a) C(T,P)/T of CeAi3 ||(for T > 8K we took the running average of 
the data points to reduce the scatter), (b) C(T,B) of CeCus.sAuo.s [ pdf| . 
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